Temperature Dependence of Standard Model CP Violation 
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We analyze the temperature dependence of CP violation effects in the Standard Model by deter- 
mining the effective action of its bosonic fields, obtained after integrating out the fermions from the 
theory and performing a covariant gradient expansion. We find non-vanishing CP violating terms 
starting at the sixth order of the expansion, albeit only in the C odd/P even sector, with coefficients 
that depend on quark masses, CKM matrix elements, temperature and the magnitude of the Higgs 
field. The CP violating effects are observed to decrease rapidly with temperature, which has impor- 
tant implications for the generation of a matter-antimatter asymmetry in the early Universe. Our 
results suggest that the cold electroweak baryogenesis scenario may be viable within the Standard 
Model, provided the electroweak transition temperature is at most of order 1 GeV. 

PACS numbers: ll.30.Er, 98.80.Bp 



Introduction. — One of the most fundamental observa- 
tions in modern cosmology is the existence of an asym- 
metry between matter and antimatter in the Universe. 
In order to generate such a state, high energy processes 
must have been efficient in breaking the baryon number, 
C (charge conjugation) and CP (combined charge conju- 
gation and parity) symmetries in an environment out of 
thermal equilibrium [H . All these ingredients are present 
in the Minimal Standard Model (SM), but the precise 
mechanism to generate the correct abundance (if one ex- 
ists) has so far eluded us 

Within the SM, CP violation originates from the com- 
plex phase of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix. One approach to include its effects in out-of- 
equilibrium field dynamics is to integrate the fermions 
out from the theory and recover CP violating bosonic 
operators in the resulting effective action — a technique 
particularly well suited for classical numerical simula- 
tions. This was first carried out in a perturbative setting 
by Shaposhnikov et al. [3], who argued that at tempera- 
tures T above the electroweak scale, CP violating opera- 
tors are suppressed by coefficients proportional to T~^^ . 
The conclusion was that the SM electroweak phase tran- 
sition, taking place at T ~ f 00 GeV, cannot constitute a 
viable mechanism for baryogenesis (see also Rcf. I^l)- 

Somewhat later, it was argued that electroweak baryo- 
genesis might still work if there was a mechanism to delay 
the transition until a temperature scale of T ~ 1-f GeV, 
a scenario typically referred to as cold electroweak baryo- 
genesis Refs. [il-fl^ subsequently employed a co- 
variant gradient expansion to determine the leading CP 
violating terms in the bosonic SM effective action at 
zero temperature; however, with contradictory results. 
In particular, while Ref. claimed the existence of a C 
even/P odd operator at the sixth order of the expansion, 
Ref. [3| reported that only the C odd/P even sector gives 
a non-vanishing CP violating contribution. Following 



these works, numerical simulations employing the zero 
temperature CP violating operator of Ref. [10| demon- 
strated that one obtains a baryon asymmetry exceeding 
the observed one by four orders of magnitude Tak- 
ing into account the known suppression of CP violation 
in the limit of very high temperatures Q , this highlights 
the importance of determining the temperature depen- 
dence of CP violating effects within the SM. 

In the present Letter, our objective is simple: To eval- 
uate the CP violating part of the SM bosonic effective 
action as a function of temperature. The resulting rate of 
temperature dependence should, in combination with the 
results of Ref. [l3| , give at least a rough estimate of the 
temperature range, in which SM cold electroweak baryo- 
genesis might be viable. In the limit of zero temperature, 
our calculation should also provide an independent check 
of the very non-trivial predictions of Ref. ll| , addressing 
the discrepancy between the results of Refs. [lol [lH. 

Bosonic effective action. — Consider a Euclidean field 
theory of chiral fermions, denoted collectively as 'ip{x) = 
{'4)'l[x),iP'b.[x)), coupled to gauge fields Vi,{x) ,Vb.{x) and 
possibly to additional scalar fields. The dynamics of 
the fermions is described by a Euclidean Lagrangian of 
the Dirac type, 'ip{x)D{x)'4>{x), where the Dirac opera- 
tor D{x) in the chiral (left-right) basis has the generic 
form 



D{x) 



+ 

■mviL{x) 



rriLRix) 



using the Feynman slash notation, f 



(1) 



The func- 



tional integral over the fermionic fields can be performed 
exactly, resulting in an additional contribution to the ef- 
fective action of the bosonic fields, 



r[^L, VR,TOLR,mRL] = -TrlogD. 



(2) 



As demonstrated in Refs. [14|, the P even and P odd 
parts of this effective action, F"*" and F~ , can be expressed 
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in the form 

r+ = -iReTr(loglf), (3) 
r- = -iiImTr(75log/<:)+rgwzw, (4) 

where the second order differential operator K is defined 
a.s K = mLRTORL - -^LL^j^L^RRmRL. The gauged 
Wess-Zumino-Witten term Fgwzw encodes the chiral 
anomaly. 

The trace-logarithm of the Dirac operator cannot 
be evaluated in a closed form for general coordinate- 
dependent background fields. Therefore, one usually re- 
sorts to a covariant gradient expansion, that is, an expan- 
sion in gauge fields and covariant derivatives. A useful 
tool for this is the method of covariant symbols 
a general method for the evaluation of traces of differen- 
tial operators of the type f{D,m), where / is an alge- 
braic function of its two arguments, a (matrix covariant) 
derivative D and a (matrix) field m{x). At nonzero tem- 
perature, the trace of this operator reads 



Tr /(-D, m) 



tr[/(D,m)l] 



(5) 



I tr[(i^oVo)V(Am)l]. 



where V,, = 



matrix 

V 



indices, and the 

ip-x,, „ip x„—iD-V 



tr" stands for a trace over 
bar denotes a conjugation. 
Finally, the inte- 
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gration measures in a d-dimensional Euclidean space- 
time have the forms = J^^"^ dxo J d'^~^x and = 

TJ2p„ J d'^"V/(27r)'^-\ where po takes the values 2n7rT 
or (2n -|- l)7rT with integer n, depending on whether the 
trace is taken over a bosonic or fermionic space. 

A detailed derivation of Eq. ([5|) will be given else- 



where 17[; here, we simply remark that at zero tem- 



perature, the derivative D as well as the field m are 
Lorentz covariant, the first term of this expression is com- 
pletely Lorentz invariant, and the second term vanishes. 
At nonzero temperature this is no longer the case. One 
should additionally note that due to the presence of ex- 
plicit free temporal derivatives in the second term, the 
individual terms in the sum are not gauge invariant. 

Application to the Standard Model. — We need not write 
down the full SM Lagrangian, but simply display its Eu- 
clidean Dirac operator in order to fix our conventions, 
following closely the notation of Ref. Upon diago- 

nalizing the mass matrices, the chiral components of the 
Dirac operator in the quark sector take the form 




„ + + \^^V \ 
U)d + 0j 



(6) 
(7) 



Here, Z^, and are the weak intermediate bo- 
son and gluon fields, respectively, and V is the CKM 
matrix. Furthermore, D^^ — + (2/3)i?p and Dd^ — 
d^i — (1/3)5^, where is the hypercharge gauge field. 
(Note that we are using a "mixed" basis in terms of 
and the hypercharge field -B^ rather than the commonly 
used Zfj^ and the photon field .) The indicated 2x2 ma- 
trix structure of the Dirac operator corresponds to weak 
isospin; Dirac, color and family indices are suppressed. 
The left-right and right-left components of the Dirac op- 
erator are identical, ttilr = torl = (</>/«) diag(m„, m^), 
where TO„,d are the (diagonal) mass matrices for u and 
d type quarks. Throughout the calculation, we use the 
unitary gauge, in which the only physical scalar field of 
the SM — the Higgs field (j) — fluctuates around its vacuum 
expectation value, v « 246 GeV. 

The above expressions for the Dirac operator com- 
pletely fix our conventions. Our fields are related to the 
canonically normalized (Euclidean) fields, denoted here 
by a tilde, via = ig/V2)W^, = [5/(2 cos^w)]^^, 
= g'Sp, and = {gs/2)\aGa^. Here, as usual, 
g, 5', 5s denote respectively the weak isospin, hypercharge 
and strong coupling constants, while is the Weinberg 
angle and \a are the Gell-Mann matrices. 

The elements of the CKM matrix are known to a good 
accuracy [3| and are fixed up to independent rephas- 
ings of the quark fields, which must leave all physical 
observables unchanged. The simplest rephasing invari- 
ant constructed from the CKM matrix is the Jarlskog 
invariant J, defined by lm{VijV^^^VMV^^) — Jeikf-ji- 
Any CP violating bosonic operator thus has to contain 
at least four fields. As shown by Smit Q, there 
is, however, no CP violation even at fourth order of 
the gradient expansion at zero temperature — a conclu- 
sion straightforwardly generalizable to nonzero tempera- 
ture. This in particular implies that the anomalous Wess- 
Zumino-Witten term need not be considered, as it only 
contributes at fourth order [ll|. Hence, the CP violating 



part of the bosonic effective action can be determined 
solely using the Dirac operators ^ and Q as well as 
Eqs. dSl) and Q. 

Finally, we note that the above exercise can naturally 
be repeated for the lepton sector of the SM, assuming the 
neutrinos to have Dirac type masses. Their contribution 
is, however, always suppressed by many orders of magni- 
tude with respect to the quarks, and we have therefore 
not considered it in our work. 

Results. — We have implemented the gradient expan- 
sion and performed the traces indicated in Eqs. and 
(j4|) using the Mathematica package Feyncalc [18|]. At 
the leading non- vanishing (sixth) order, the CP violating 
part of the effective action can be written in the form 



Fcp-o. 



dd 



(8) 



-'-N.JGfkcp dxojd^x{^ (O0 + O1+O2), 
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where iVc = 3 is the number of colors, Gp = l/(\/2w^) 
the Fermi couphng, and the zero temperature effective 
couphng Kcp reads 



KCP 



4(p')'^7;;CT^«3.lxlo^ 



x(TO^-TO^)(m^-77i^)(m2 -mg) 



(9) 



The index n in the functions ©„ counts the number of 
Z or tp fields where (p^ = Each of these three 

terms can be subsequently divided into P even and P odd 
parts, On = +0~ . These functions contain not only 
Lorentz invariant operators, but also ones containing one 
or several vectors = S^q, specifying the rest frame of 
the thermal bath. The non-invariant terms must natu- 
rally vanish in the limit of zero temperature, which we 
have indeed verified. 

The Lorentz invariant operators in read explicitly 



3 

Cl 

3 

2ci 



5C2 , 



-—W^W^W~ w~ 

+ ^2^ 



C.C., (10) 



-c^{W+ ■ W-)W+W~, 

+cj{W+ ■ W-)W+W-, + c^W+W+W-W-, 

~c,2{w+ ■ w~)w+w^^ - c,2W+w+w~w;;^ 

+ci:iW;W+W+W-^] - C.C., (11) 



X [cs{W+fW;W- ~ cs{W-fW+W+ 
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{z ■ p) [cg{w+ -w-y- 2cq{w+Y{w-) 



X [cw[W+fW;W- + cw{W-fW+W+ 

-2cn(w^+ • w~){w+w;; + w+w;:)] , (12) 

while we find that 



On 



or = = 0. 



(13) 



Here, "c.c." stands for complex conjugation, acting on 
the fields as -J> -W^, Z -> -Z, and ^ ^p. In ad- 
dition, we have denoted the hypercharge covariant deriva- 
tives of by W^j^ = {d^ ± B^)Wj^. For the numerical 
values of the quark masses and the Jarlskog invariant (see 
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FIG. 1: The coefficients ci-cis plotted as functions of the 
effective temperature Tcfr = vT/<j). The bold fines correspond 
to the smallest and largest of the d that approach one at zero 
temperature, i.e. Ci and cio. 



Ref. [ig|), we have used m„ = 2.5 MeV, = 5 MeV, 
mc = 1.27 GeV, = 100 MeV, = 172 GeV, 
nib = 4.2 GeV as well as J = 2.9 x 10"^. 

The coefficients Ci depend on the quark masses to/ as 
well as the temperature T and the Higgs field 4'{x), which 
appear in the particular combination Teg = vT/cjj. In the 
zero temperature limit, ci-cn approach unity while Ci2 



Hi 



ci3 tend to zero, reducing our result to that of Ref. 
and thus independently verifying its conclusions. 

In Fig. [U we display the behavior of the coefficients c; 
as functions of Toff, obtained through numerical evalua- 
tion of their defining one-loop, multi scale sum-integrals. 
At low temperatures, the Ci-Cn evolve slowly until at 
Teff ~ 10-30 MeV they begin to fall rapidly. The C12-C13 
on the other hand first exhibit a fast increase, reaching 
their maximum around 20-30 MeV, after which they, too, 
begin to decrease. The largest coefficient at all temper- 
atures is cio, which reaches its maximal value of 1.3 at 
Toff « 18 MeV. In addition, we note that at Toff > 10 
GeV, all Ci's are at most of order 10-^"*^, consistent with 
the estimates of Shaposhnikov et al. Q . 

Fig. [21 on the other hand, demonstrates the depen- 
dence of our results on the value of the Higgs field (j). 
As is evident from the functional form of Toff = vT/(j), 
larger values of (j) ameliorate the thermal suppression of 
CP violating effects. This observation has important im- 
plications for cold electroweak baryogenesis simulations, 
highlighting the necessity of determining the distribution 
of the Higgs field during the cold spinodal transition . 

Finally, we want to stress that in the present work we 
have completely neglected the effects of the strong inter- 
action, an issue of increasing severity as one approaches 
the deconfinement transition (see also the discussion in 
Ref. [l^). For an extended discussion of this issue as 
well as of the details of our computation (including the 
Lorentz breaking operators present at nonzero tempera- 
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FIG. 2: The dependence of the d on the Higgs field 0, plotted 
for three different temperatures, T = 1/3 GeV, 1 GeV and 
3 GeV. In each case, the grey band is spanned by ci and cio- 



ture), we refer the reader to Ref. [171 ]. 

Conclusions and outlook. — In this Letter, we have de- 
termined the leading CP violating operators in the co- 
variant gradient expansion of the bosonic effective ac- 
tion of the Standard Model, displayed in Eqs. ((St- P^ 
above. In the zero temperature limit, we independently 
confirmed the result of Ref. 



11| . While this agreement 



itself does not resolve the conflict between the results of 
Refs. 0, 
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it suggests that the source of the discrep- 
ancy must lie either in the details of the computation 
performed in Ref. 10], or in a subtlety in separating T~ 
into an anomalous and a regular part, cf. Eq. 

The gradient expansion we employ involves setting the 
external momenta to zero. One should bear in mind that 
this may limit the range of applicability of the result. 
In agreement with Refs. [13, U], we expect the gradi- 
ent expansion to be valid up to external momenta of the 
order of nic- This is based on the observation that the 
result is insensitive to the three smallest quark masses. 
In Ref. [l3] , we address the issue of convergence by cal- 
culating the next order in the expansion. 

Nonzero temperature was observed to affect the values 
of the couplings of the operators already present at zero 
temperature, as well as to give rise to altogether new 
operators. Still, similarly to the zero temperature limit, 
all CP violating operators turned out to be P even/C 
odd. Another important feature of our result was that 
the temperature entered only through the combination 
Tcff = vT/cf), involving the local Higgs field </> and its 
vacuum expectation value v. 

Direct numerical simulations using the SM effective ac- 
tion with a P-odd operator (of the O^-type in our classi- 
fication) suggest that the cold electroweak baryogenesis 
scenario is able to produce the observed baryon asym- 
metry of the Universe with Ci ~ 10~^ [13, [2l|. During 
the phase transition, the fields are strongly out of equi- 
librium, and thus one would ideally like to derive the CP 



violating operators in such an environment. As this task 
is, however, technically extremely demanding (unless one 
includes the fermions themselves in the dynamics, which 
has become viable only recently [l^), our calculation as- 
sumed thermal equilibrium at an effective temperature T. 
As a rough estimate of the temperature range, in which 
cold electroweak baryogenesis might be viable, we ob- 
serve that for (j) — V, the largest Ci coefficient ciq reaches 
the value 10^^ at T ~ 1 GeV. It should be noted, how- 
ever, that the effective action of Eqs. (j8|)~ ([T2|) consists of 
a variety of operators with coefficients of different mag- 
nitude, which may either enhance each others' effects or 
lead to surprising cancelations. 

The result we have obtained should be contrasted with 
the energy scales typically associated with cold elec- 
troweak baryogenesis. If the relevant momentum scale 
of the dynamics is identified with the inverse time scale 
of the electroweak transition, lattice simulations using 
a hypothetical order four operator indicate an effective 
temperature of T ~ Tq^ ~ 5-10 GeV, where tq is the 
maximal time scale of a transition sufficiently out of 
equilibrium to generate a large enough baryon asymme- 
try [23|. These numbers are somewhat, but not dramati- 
cally higher than the approximative 1 GeV we obtained. 
To further quantify the result, classical numerical simu- 
lations of the bosonic effective action, including all of the 
operators we derived, are clearly required. 

An important technical note in this context is that, 
since in the bosonic simulations the asymmetry is quanti- 
fied via the Chern-Simons number rather than the baryon 
number itself, the dynamics have to be explicitly CP but 
also P violating to generate an asymmetry. Given that 
the CP violating terms found here conserve P, one must 
in addition include the C odd/P odd but CP even opera- 
tors. These are inherent in the fermionic sector of the SM 
and appear at lower orders in the derivative expansion of 
the bosonic effective action than the CP violating op- 
erators. We will further address this important point in 
Ref. [l7|. Beyond this, improving the accuracy of the pre- 
dictions would require not integrating out the fermions 
but rather incorporating them into the (quantum) sim- 
ulations. We plan to return to these questions in future 
work. 
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